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Abstract. The Time-dependent Hartree-Fock (TDHF) theory is applied to the study of heavy 
ion collisions at energies around the Coulomb barrier. The competition between fusion and 
nucleon transfer mechanisms is investigated. For intermediate mass systems such as 16 O+ 20S Pb, 
proton transfer favors fusion by reducing the Coulomb repulsion. A comparison with sub-barrier 
transfer experimental data shows that pairing correlations are playing an important role in 
. enhancing proton pair transfer. For heavier and more symmetric systems, a fusion hindrance 

is observed due to the dominance of the quasi- fission process. Typical quasi-fission time of 
few zeptoseconds are obtained. Actinide collisions are also investigated both within the TDHF 
approach and with the Ballian-Veneroni prescription for fluctuation and correlation of one-body 
observables. The possible formation of new heavy neutron-rich nuclei in actinide collisions is 
discussed. 
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1. Introduction 

Microscopic description of nuclear reactions at low energy where quantum effects play a 
significant role is an important challenge of nuclear physics. The interplay between nuclear 
structure and reaction mechanisms is crucial at energies around the Coulomb barrier generated 
by the competition between Coulomb and nuclear interactions. For instance, this has been well 
established for fusion reactions (for a review, see, e.g., Ref. [1]). 

A good understanding of these effects is crucial for the quest of superheavy elements (SHE) [2] . 
Indeed, different structures of the collision partners may change SHE production cross sections 
by several orders of magnitudes. For instance, collisions of quasi-symmetric heavy systems with 
typical charge products Z1Z2 > 1600 exhibit a fusion hindrance when compared to systematics 
and models established for lighter systems [3]. The dynamical shape evolution of the dinuclear 
system and the transfer between the fragments are expected to play a significant role on this 
fusion hindrance. 

In the case of the collision of actinide nuclei, the Coulomb repulsion is so strong that fusion 
cannot occur. However, a dinuclear system may survive up to few zeptoseconds, forming the 
heaviest nuclear systems available on Earth. Multi-nucleon transfer between actinides may 
provide an alternate way to produce more neutron rich heavy and super-heavy elements. The 



complex dynamics of such systems with up to ~ 500 nucleons in interaction is very challenging 
and has been investigated with different theoretical approaches recently [4, 5, 6, 7, 8, 9, 10, 11, 12]. 

The time-dependent Hartree-Fock (TDHF) theory [13] presents the advantage of treating 
both structure and reaction dynamics within the same formalism. The TDHF theory provides 
a self-consistent mean-field description of the many-body dynamics. Early TDHF calculations 
in nuclear physics used various symmetries and simplified Skyrme interactions [14] to reduce 
computational times [15, 16, 17]. Recent increase of computational power allowed realistic TDHF 
calculations in 3 dimensions [18, 19, 20, 21, 22] with full Skyrme energy density functional (EDF) 
including spin-orbit terms [14, 23]. Several examples of recent applications of TDHF to nuclear 
dynamics were presented at the NN2012 conference [24, 25, 26, 27]. 

In Section 2, we recall the TDHF formalism and provide some details on the calculations. 
In Section 3, we study the interplay between fusion and transfer in the 16 O+ 208 Pb system. 
Section 3 is dedicated to the description of sub-barrier transfer in this reaction. In Section 5, we 
investigate the fusion hindrance of heavy quasi-symmetric systems and the quasi-fission process. 
Finally, we study the collision dynamics of actinide nuclei in Section 6 before to conclude. 

2. The Time-Dependent Hartree-Fock approach 

In nuclear physics, the TDHF theory is applied with a Skyrme EDF modeling nuclear interactions 
between nucleons [14]. The EDF is the only phenomenological ingredient of the model, as it has 
been adjusted on nuclear structure properties like infinite nuclear matter and radii and masses 
of few doubly magic nuclei [23] . The main approximation of the theory is to constrain the many- 
body wave function to be an antisymmetrized independent particle state at any time. The latter 
ensures an exact treatment of the Pauli principle during time evolution. Though TDHF does 
not include two-body collision term, it is expected to treat correctly one-body dissipation which 
is known to drive low energy reaction mechanisms as Pauli blocking prevents nucleon-nucleon 
collisions. 

The main advantage of TDHF is that it treats static properties and dynamics of nuclei 
within the same formalism and the same EDF. The initial state is obtained through static HF 
calculations which give a good approximation of nuclear binding energies and deformations [28, 
29]. Another important advantage of TDHF for near-barrier reaction studies is that it contains 
all types of couplings between the relative motion and internal degrees of freedom. However, 
TDHF gives only classical trajectories for the time-evolution and expectation values of one-body 
observables. In particular, it does not include tunneling of the many-body wave function and 
may underestimate fluctuations of one-body observables [30, 31]. The latter fluctuations can 
be estimated in the time-dependent RPA limit using a prescription proposed by Balian and 
Vcncroni [32]. 

Inclusion of pairing correlations responsible for superfluidity in nuclei have been done recently 
to study vibrations in nuclei [33, 34, 35, 36]. However, realistic applications to heavy ion collisions 
are not yet achieved and are beyond the scope of this work. 

2.1. Formalism 

The TDHF equation can be written as a Liouville-Von Neumann equation 



where p is the one body density matrix associated to the total independent particle state with 
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The sum runs over all occupied single particle wave functions <pi and r, s and q denote 
the position, spin and isospin of the nucleon, respectively. The Hartree-Fock single particle 
Hamiltonian h[p] is related to the EDF, noted E[p], by its first derivative 
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2.2. Practical aspects 

A TDHF calculation of two colliding nuclei is performed assuming that the two collision partners 
are initially at a distance Dq in their HF ground state. This distance has to be large enough 
to account for Coulomb excitation in the entrance channel (polarization, vibration, rotation...). 
Typical initial distances are within the range Dq ~ 30 — 50 fm. We assume that before the 
initial time, the nuclei followed a Rutherford trajectory. This assumption determines the initial 
velocities in the center of mass frame. The collision partners are put into motion using a Galilean 
boost [37] which is applied on each nucleus at the first iteration. 

We use the tdhf3d code built by P. Bonche and coworkers with the SLy4<i Skyrme 
parametrization [18]. This code has a plane of symmetry (the collision plane). It uses the 
Skyrme energy functional expressed in Eq. (A. 2) of Ref. [38] where the tensor coupling between 
spin and gradient has been neglected. The lattice spacing is Ax = 0.8 fm and the time step is 
At = 1.5 x 10~ 24 s. More details on modern TDHF calculations to nuclear dynamics can be 
found in Refs. [39, 40]. 

3. Fusion with medium mass systems 

A natural application of the TDHF approach in nuclear physics is to study fusion reactions at 
and above the barrier. Fusion occurs by transferring relative motion into internal excitation via 
one-body mechanisms well treated by the TDHF approach. As a result, early TDHF codes have 
been successfully applied to describe above-barrier fusion reaction in light systems [41]. In this 
section, we investigate the fusion process in the 16 O+ 208 Pb system. 

A reference nucleus-nucleus potential could be obtained in the frozen approximation with 
HF (or HFB) densities [42], where the collision partners are assumed to keep their ground-state 
density during the approach. The frozen potential can be computed with the same Skyrme 
EDF as in the TDHF calculations by translating the nuclei in their HF state. A comparison 
between TDHF and frozen fusion barriers allows to identify the role of dynamical effects, which 
are included in TDHF but absent in the frozen approach. 

Writing the HF energy E[p] as an integral of an energy density H[p(r)], i.e., 



where R is the distance between the centers of mass of the nuclei, and p\p are the densities of 
their HF ground-state. Eq. (5) neglects the Pauli principle between the nucleons of one nucleus 
and the nucleons of the other one. As a result, it is valid only for small overlaps between the 
two collision partners. However, for light and intermediate mass systems, the barrier radius is 
large enough to justify this approximation. 

The HF-frozen potential for the 16 O+ 208 Pb system is shown in Fig. 1. The potential obtained 
from the Wong formula [43] is also presented for comparison. The height of the barrier obtained 
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Figure 1. Nuclear (dashed blue line) and Coulomb (dotted red line) contributions of 
the nucleus-nucleus potential (solid line) obtained with the HF-frozen approximation in the 
16 O+ 208 Pb system. The Wong potential [43] is shown in green dot-dashed line. The latter is 
obtained with a potential depth Vq = 70 MeV, a potential diffuseness a = 0.48 fm, and nuclear 
radii Rj = 1.25A 1 / 3 fm. 



with the HF-frozen approximation is y^ rozen ~ jqq MeV at Rg° zen ~ 11.8 fm. Note that the 
same barrier height has been obtained in Refs. [44, 45] with a similar implementation of the HF- 
frozen approximation. This value is close to the barrier obtained with the Wong formula [43], 

yWong _ ?5 g MeV> whne it ig 1 smaller than the Bags barrier [4 6 ] ; yBass ^ J7 q Me y 

at Rg ass ~ 11.4 fm. A comparison with the experimental barrier distribution shown in Fig. 2 
indicates that all these barriers overestimate the experimental value V^ xp ' ~ 74.5 MeV obtained 
from the centroid of the barrier distribution. 

We now investigate the possible role of dynamical effects on the fusion barrier with TDHF 
calculations. The TDHF fusion barrier is identified as the capture threshold for central collisions, 
above which a compound system is formed and below which two fragments are emitted. Due 
to the finite time of the TDHF evolutions, one has to define a maximum computational time 
(~ 10 3 fm/c in the present case) above which the final configuration (i.e., one compound system 
or two fragments) is assumed to be reached. 

Figure 3 shows the evolution of the relative distance between fragment centers of mass in 
central 16 O+ 208 Pb collisions at different energies around the barrier. We identify two sets of 
trajectories associated to fusion for E > 74.45 MeV and to re-separation of the fragments 
for E < 74.44 MeV. The resulting fusion barrier predicted by these TDHF calculations is, 
then, Vb DHF = 74.445 ± 0.005 MeV. As a result, the dynamical effects included in TDHF 
calculations lower the barrier by ~ 1.5 MeV for this system as compared to the HF-frozen 
approximation. Other methods based on a macroscopic reduction of the mean-field dynamics, 
namely the dissipative-dynamics TDHF [44] and the density-constrained TDHF [48], also find 
similar results with an energy dependence to the barrier heights ranging from 74.5 at low energies 
to 76 MeV at higher energies where the frozen approach is more valid. 

The density evolutions close to the fusion threshold, i.e., at E cm . = 74.44 and 74.45 MeV, 
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Figure 2. Experimental fusion barrier distribution of the 16 O+ 208 Pb system from Ref. [47]. 
The HF-frozen, Bass and TDHF barriers are indicated by arrows. 




Figure 4. (top) Density evolution for the reaction 16 O+ 208 Pb corresponding to a head-on 
collision at a center of mass energy E c _ m _ = 74.44 MeV (just below the fusion barrier). The red 
surfaces correspond to an iso-density at half the saturation density (po/2 = 0.08 fm~ 3 ). Each 
figure is separated by a time step of 135 fm/c. Time runs from left to right, (bottom) Same at 
E c .m. = 74.45 MeV, i.e., just above the fusion threshold. 

are plotted in Fig. 4. In both "di- nuclear" system is formed during a relatively long 

time (~ 500 fm/c) before either re-separation or fusion. Nucleon transfer is expected to occur 
in this di-nuclear system. An evidence of this transfer is given by integrating the proton and 
nucleon densities at the final time when the two fragments re-separate at E c . m . = 74.44 MeV. As a 
result, the outgoing channel is, in average, 14 C+ 210 Po, showing the importance of proton transfer 
from the light to the heavy partner, in good agreement with experiment [49, 50]. This proton 
transfer effectively lowers the barrier by decreasing Z1Z2 and, then, the Coulomb repulsion. As a 
result, the coupling to transfer channels increase the fusion probability for this system. Transfer 
reactions in the 16 O+ 208 Pb system are discussed in more details in the next section. Note that 
low- lying collective vibrations are also known to affect the fusion barrier distribution [47]. 

4. Transfer reactions 

Transfer reactions in heavy-ion collisions have been investigated within the TDHF framework 
recently [51, 52, 53, 39, 54, 55]. In fact, the previous study showed the importance of the 
interplay between fusion and transfer reactions in the 16 O+ 208 Pb system around the barrier. 
For instance, just below the barrier, the light fragment in the exit channel (see top of Fig. 4) is, 
in average, a 14 C nucleus, indicating that the transfer of two protons is the dominant channel 
at this energy. We now discuss this sub-barrier transfer reaction in more details. 

In the previous section, a signature for transfer mechanism was obtained from the change of 
average numbers of nucleons in the fragments. Using the quantum nature of the TDHF theory, 
it is possible to compute the probability for a specific outgoing channel, i.e., the probability to 
find a given number of neutrons and protons in the fragments [56, 54]. This is possible, e.g., 
thanks to a particle number projection technique applied on the outgoing fragments [54]. 

Let us consider a final time when the two fragments are located on both side of the x = 
plane and define a particle number projector onto TV" protons or neutrons in the x > region 
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where 
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counts the number of particles in the x > region (B(x) = 1 if x > and elsewhere). Isospin 
is omitted to simplify the notation. The projector defined in Eq. (6) can be used to compute 
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Figure 5. Neutron (circles) and proton (squares) number probability distributions of the lightest 
fragment in the exit channel of a head-on 16 O+ 208 Pb collision at E c . m , = 74.44 MeV (left) and 
65 MeV (right). Adapted from Ref. [54]. 



the probability to find N nucleons in x > in the final state \(j>), 

P R (N)\<f>) 2 = (<t>\P R (N)\<t>) = ±- / 2 d0 e-Wtfltoie)), (8) 
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where \4> R {9)) = e R \4>). Note that \4>r{9)) is an independent particle state and, then, the last 
term in Eq. (8) is the determinant of the matrix of the occupied single particle state overlaps [54]: 



b R (6)) = det(F) (9) 
with 

F lj = Y: /dr ^*(rV|(r)e^ e W (10) 

s 

The integral in Eq. (8) is discretized using 9 n = 2irn/M with the integer n = 1 • • • M. Choosing 
M = 300 ensures numerical convergence for the 16 O+ 208 Pb system. Fig. 5 shows the resulting 
transfer probabilities at E c _ m4 = 74.44 MeV (left) and at E c .m. = 65 MeV (right). In agreement 
with the results presented in section 3, the most probable channel is Z = 6 and N = 8 at the 
barrier. However, lowering the energy reduces the transfer probabilities and the main channel 
is Z = N = 8 well below the barrier, corresponding to inelastic and elastic channels. 

To compare with experimental data on the 16 O+ 208 Pb reaction, we plot in Fig. 6 the transfer 
probabilities from TDHF as a function of the distance of closest approach R m in between the 
collision partners [57], assuming a Rutherford trajectory [58]: 

R min = Z x Z 2 e 2 \\ + cosec(0 c . m ./2)]/2£ c . m . (11) 

where 9 c _ m _ is the center of mass scattering angle. Recent data from Ref. [50] are shown in 
Fig. 6 for sub-barrier one- and two-proton transfer channels. We see that TDHF overestimates 
the one-proton transfer probabilities and underestimates the strength of the two-proton transfer 
channel. This discrepancy is interpreted as an effect of pairing interactions [54, 50]. Indeed, 
paired nucleons are expected to be transferred as a pair, increasing (resp. decreasing) the 
two-nucleon (single-nucleon) transfer probability. For R m in > 13 fm, however, the TDHF 
calculations reproduces reasonably well the sum of one and two-proton transfer channels. For 




Figure 6. Proton number probability as function of the distance of closest approach in the small 
outgoing fragment of the 16 O+ 208 Pb reaction. TDHF results are shown with lines. Experimental 
data (open symbols) are taken from Ref. [50]. 



Rmin < 13 fm, sub-barrier fusion, not included in the TDHF calculations, reduces transfer 
probabilities [54, 50]. 

These studies emphasize the role of pairing interactions in heavy-ion collisions. The recent 
inclusion of pairing interactions in 3-dimensional microscopic codes [34, 35, 36] gives hope in our 
ability to describe such data more realistically in a near future. 

5. Fusion and quasi-fission in heavy systems 

To a reasonably good approximation, the fusion barrier for light and intermediate mass systems 
is determined by the frozen barrier. This approximation fails, however, for heavy systems with 
typical charge products Z1Z2 > 1600 which are known to exhibit fusion hindrance [3]. Above 
this threshold, an extra-push energy is usually needed for the system to fuse [59]. In fact, at 
the energy of the frozen barrier, heavy systems are more likely to encounter quasi-fission, i.e., a 
re-separation in two fragments after a possible mass exchange. 

We now investigate the reaction mechanism in heavy systems with possible fusion hindrance. 
We first illustrate the fusion hindrance with TDHF calculations of fusion thresholds. Then, we 
investigate the quasi-fission process. 

5.1. TDHF calculations of fusion hindrance 

Let us first consider the 90 Zr+ 124 Sn system which has a charge product Z1Z2 = 2000, and, then, 
is expected to exhibit a fusion hindrance. Indeed, the proximity model [60] predicts a barrier for 
this system V prox ' ~ 215 MeV, while TDHF calculations predicts that the system encounters a 
fast re-separation at this energy, as shown in Fig. 7 [61]. In the same figure, we observe a long 
contact time at E c , m . = 240 MeV, which is interpreted as a capture trajectory leading to fusion. 
The additional energy needed to fuse is then ~ 25 MeV, which is higher than the extra-push 
model [59] prediction E x ~P ush ~ 14.8 MeV. 

TDHF calculations of fusion hindrance have been performed for several systems, as shown in 
Fig. 8, where the TDHF fusion thresholds are compared with the interaction barriers predicted 
by the proximity model [60] and with the extra-push model [59]. We observe an increase of 




Figure 7. Distance between the centers of mass of the fragments as a function of time in 
head-on 90 Zr+ 124 Sn collisions for different center of mass energies. 
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Figure 8. TDHF fusion thresholds for several heavy systems are compared with the proximity 
barrier [60] and with results from the extra-push model [59]. 



the fusion threshold with TDHF as compared to the proximity model which assumes frozen 
reactants. This indicates that dynamical effects are playing an important role in the reaction 
by hindering fusion. The order of magnitude of the additional energy needed to fuse is similar 
to the one predicted with the extra-push model. 

5.2. The quasi-fission process 

The quasi-fission mechanism is a fast re-separation of the fragments, with usually a partial 
equilibration of their mass. Typical quasi-fission times are of the order of few zepto- 
seconds [62, 63, 40]. Quasi-fission becomes dominant in heavy systems and is mostly responsible 
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Figure 9. Density profile in the Zr+ Sn head-on collision at E c/m . = 235 MeV. From 
Ref. [61]. 
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Figure 10. Snapshots of the TDHF isodensity at half the saturation density in the Ca+ U 
system for different initial orientations and E c . m _. 



for the fusion hindrance discussed in the previous section. As an example, Fig. 9 shows the 
density profiles for the 90 Zr+ 124 Sn head-on collision at E c . m . = 235 MeV. We observe that 
the two fragments are in contact during ~ 5 zs. A dinuclear system is then formed before 
re-separation in two fission- like fragments. 

Extensive calculations are ongoing on the 40 Ca+ 238 U system to compare with recent 
measurements performed at the Australian National University [64]. Examples of density 
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Figure 11. TDHF calculations of the mass of the fragments (top) and of the quasi- fission time 
(bottom) in 40 Ca+ 238 U central collisions as a function of the center of mass energy (divided by 
the proximity barrier [60]). For quasi-fission times larger than 23 zs, only a lower limit is given. 
Two different orientations of the 238 U are considered (see inset). 



evolutions are shown in Fig. 10 for two different initial conditions. In both cases, a quasi- 
fission process is obtained. In particular, we observe an important multi-nucleon transfer from 
the heavy fragment toward the light one. 

Fig. 11 presents final fragment masses (top) and quasi-fission times (bottom) for two different 
orientations of the 238 U. Comparing these two figures, we observe that the mass equilibration 
(i.e., the formation of two fragments with symmetric masses) is not complete and varies with 
the life-time of the dinuclear system, i.e., the longer the contact, the larger the mass transfer. 
We also see that all the calculations with the 238 U deformation axis aligned with the collision 
axis lead to a quasi-fission with partial mass equilibration and quasi-fission times smaller than 
10 zs. Shell effects may affect the final outcome of the reaction by favouring the production 



Figure 12. Snapshots of the isodensity at half the saturation density in 2cS8 u+ U central 
collisions at E c , m , = 900 MeV from TDHF calculations. Snapshots are given at times t = 0, 1.5, 
2.7, and 4.2 zs from top to bottom. From Ref. [11]. 



of fragments in the Pb region. In particular, this orientation never leads to fusion, while 
the other orientation produces long contact times above the barrier which may be associated to 
fusion. Calculations of non-central 40 Ca+ 238 U collisions are ongoing in order to compare with 
experimental data. 

6. Actinide collisions 

Collisions of actinides form, during few zs, the heaviest nuclear systems available on Earth. In 
one hand, such systems are interesting to study the stability of the QED vacuum under strong 
electric fields [65, 66, 11]. In the other hand, they might be used to form neutron-rich heavy 
and super- heavy elements via multi-nucleon transfer reactions [6, 7, 12]. 

Actinide collisions have been studied with the TDHF approach [67, 11, 12, 68]. Fig. 12 
shows density evolutions in 23 8tj_|_238tj cen tral collisions at E c . m . = 900 MeV for different initial 
conditions. We see that the initial orientation of the nuclei plays a crucial role on the reaction 
mechanism [11], with the production of a third fragment in the tip-tip collision (left), or net mass 
transfer in the tip-side configuration (middle). In the latter case, ~ 6 protons and ~ 11 neutrons, 
in average, are transferred from the right to the left nucleus, corresponding to the formation of 
a 255 Cf primary fragment. 

Similar calculations have been performed on the 232 Th+ 250 Cf system [12]. An example of 
density evolution is shown in the right panel of Fig. 13. In this case, we observe a net transfer 
of nucleons from the tip of the 232 Th to the side of the 250 Cf, corresponding to an inverse 
quasi-fission process, i.e., the exit channel is more mass-asymmetric than the entrance channel. 
Indeed, in this case, a 265 Lr fragment is formed in the exit channel. It is worth mentioning 
that these calculations predict inverse quasi-fission for this specific orientation only, i.e., when 
the tip of the lighter actinide is in contact with the side of the heavier one. Indeed, the other 
orientations induce "standard" quasi-fission [12]. Note also that another inverse quasi-fission 
mechanism is predicted in actinide collisions due to shell effects in the 208 Pb region [69, 6, 7]. 

In the previous example, the 265 Lr heavy fragment indicates the average N and Z of a 
distribution. The fluctuations and correlations of these distributions have been computed with 
the Balian-Veneroni prescription [32] using the tdhf3d code [70, 68]. Fig. 13(left) shows the 
resulting probabilities assuming Gaussian distributions of the form 
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Figure 13. (right) Snapshots of the isodensity at half the saturation density in 232 Xh+ 250 Cf 
central collisions at E c . m . = 916 MeV. (left) Gaussian distributions of N and Z with widths and 
correlations computed with the BV prescription (linear color scale). The solid line represents 
the predicted j3— stability line. 



where <jn,z are the standard deviations and < p < 1 quantifies the correlations between the 
N and Z distributions [40]. We see that many new f3— stable and neutron-rich heavy nuclei 
could be produced if their excitation energy is low enough to allow their survival against fission. 
This inverse quasi-fission process needs further investigations, in particular to predict realistic 
production cross-sections. 

7. Conclusions and perspectives 

The TDHF theory has been applied to the study of heavy ion collisions at energies around the 
Coulomb barrier. Its time-dependent nature allows to investigate dynamical effects responsible 
for the modification of the fusion thresholds. In particular, the coupling between the relative 
motion and proton transfer in 16 0+ 208 decreases the barrier by reducing the Coulomb repulsion. 

Particle transfer probabilities are predicted using a particle-number projection technique. A 
comparison with experimental data shows the importance of pairing correlations on transfer. The 
latter could be studied with new time-dependent Hartree-Fock-Bogoliubov codes. In particular, 
one could answer the question on the origin of these pairing correlations: Are they present in 
the ground-states or are they generated dynamically during the collision? 

For heavy and (quasi-) symmetric systems, a fusion hindrance is observed due to the quasi- 
fission process. The possibility to study the quasi-fission mechanism with a fully microscopic 
quantum approach such as the TDHF theory is promising. It will help to understand the strong 
fusion hindrance in quasi-symmetric heavy systems and may provide a guidance for new fusion 
experiments with exotic beams. 

Actinide collisions have been investigated both within the TDHF approach and with the 
Ballian-Veneroni prescription. A new inverse quasi-fission mechanism associated to specific 
orientations was found. This mechanism might help to produce new neutron-rich heavy nuclei. 
A systematic investigation of this effect is mandatory to help the design of future experimental 



equipments dedicated to the study of fragments produced in actinide collisions. 



Acknowledgements 

D. Hinde, M. Dasgupta and M. Evers are warmly thanked for useful discussions. A. W. is 
grateful to CEA/Saclay, IRFU/SPhN, where part of this work has been performed. The 
TDHF calculations were performed on the NCI National Facility in Canberra, supported 
by the Commonwealth Government. Support from ARC Discovery grants DP06644077 and 
DP110102858 is acknowledged. 



References 

M. Dasgupta, D.J. Hinde, N. Rowley, A.M. Stefanini, Ann. Rev. Nucl. Part. Sci. 48, 401 (1998) 
S. Hofmann, G. Miinzenberg, Rev. Mod. Phys. 72, 733 (2000) 

H. Gaggeler, T. Sikkeland, G. Wirth, W. Briichle, W. Bogl, G. Franz, G. Herrmann, J.V. Kratz, M. Schadel, 
K. Siimmerer et al., Z. Phys. A 316, 291 (1984) 

G.G. Adamian, N.V. Antonenko, A.S. Zubov, Phys. Rev. C 71, 034603 (2005) 
Z.Q. Feng, G.M. Jin, J.Q. Li, Phys. Rev. C 80, 067601 (2009) 

V.I. Zagrebaev, Y.T. Oganessian, M.G. Itkis, W. Greiner, Phys. Rev. C 73, 031602 (2006) 
V. Zagrebaev, A. Karpov, W. Greiner, in NN2012 Proceedings, to appear in Journal of Physics: Conference 

Series (JPCS) (2012), arXiv: 1207.5700 
T. Maruyama, A. Bonasera, M. Papa, S. Chiba, Eur. Phys. J. A 14, 191 (2002) 
J. Tian, X. Wu, K. Zhao, Y. Zhang, Z. Li, Phys. Rev. C 77, 064603 (2008) 
K. Zhao, X. Wu, Z. Li, Phys. Rev. C 80, 054607 (2009) 

C. Golabek, C. Simenel, Phys. Rev. Lett. 103, 042701 (2009) 

D. J. Kedziora, C. Simenel, Phys. Rev. C 81, 044613 (2010) 
P.A.M. Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930) 
T. Skyrme, Phil. Mag. 1, 1043 (1956) 

Y. Engel, D. Brink, K. Goeke, S. Krieger, D. Vautherin, Nucl. Phys. A 249, 215 (1975) 
P. Bonche, S. Koonin, J.W. Negele, Phys. Rev. C 13, 1226 (1976) 
J.W. Negele, Rev. Mod. Phys. 54, 913 (1982) 
K.H. Kim, T. Otsuka, P. Bonche, J. Phys. G 23, 1267 (1997) 

J.A. Maruhn, P.G. Reinhard, P.D. Stevenson, J.R. Stone, M.R. Strayer, Phys. Rev. C 71, 064328 (2005) 
T. Nakatsukasa, K. Yabana, Phys. Rev. C 71, 024301 (2005) 

A. S. Umar, V.E. Oberacker, Phys. Rev. C 71, 034314 (2005) 
F. Sebille, S. Figerou, V. de la Mota, Nucl. Phys. A 822, 51 (2009) 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. A 635, 231 (1998) 
N. Itagaki et al., in NN2012 Proceedings, to appear in Journal of Physics: Conference Series (JPCS) (2012) 
Y. Iwata, S. Heinz, in NN2012 Proceedings, to appear in Journal of Physics: Conference Series (JPCS) 

(2012), arXiv: 1208. 6215 

V.E. Oberacker, A.S. Umar, R. Keser, in NN2012 Proceedings, to appear in Journal of Physics: Conference 

Series (JPCS) (2012), arXiv: 1208 . 0263 
V. de la Mota et al, in NN2012 Proceedings, to appear in Journal of Physics: Conference Series (JPCS) 
(2012) 

D. Vautherin, D.M. Brink, Phys. Rev. C 5, 626 (1972) 
K. Bennaceur, P. Bonche, J. Meyer, C. R. Physique 4, 555 (2003) 

C. H. Dasso, T. D0ssing, H.C. Pauli, Z. Phys. A 289, 395 (1979) 
R. Balian, M. Veneroni, Phys. Rev. Lett. 47, 1353 (1981) 
R. Balian, M. Veneroni, Phys. Lett. B 136, 301 (1984) 

B. Avez, C. Simenel, P. Chomaz, Phys. Rev. C 78, 044318 (2008) 
M. Assie, D. Lacroix, Phys. Rev. Lett. 102, 202501 (2009) 

S. Ebata, T. Nakatsukasa, T. Inakura, K. Yoshida, Y. Hashimoto, K. Yabana, Phys. Rev. C 82, 034306 
(2010) 

I. Stetcu, A. Bulgac, P. Magierski, K.J. Roche, Phys. Rev. C 84, 051309 (2011) 

D. Thouless, J. Valatin, Nucl. Phys. 31, 211 (1962) 
P. Bonche, H. Flocard, P. Heenen, Nucl. Phys. A 467, 115 (1987) 

C. Simenel, D. Lacroix, B. Avez, Quantum Many-Body Dynamics: Applications to Nuclear Reactions (VDM 
Verlag, Sarrebruck, Germany, 2010), arXiv: 0806. 2714 

[40] C. Simenel, review article to appear in Eur. Phys. J. A. 



[41] P. Bonche, B. Grammaticos, S. Koonin, Phys. Rev. C 17, 1700 (1978) 

[42] V.Y. Denisov, W. Norenberg, Eur. Phys. J. A 15, 375 (2002) 

[43] C.Y. Wong, Phys. Rev. Lett. 31, 766 (1973) 

[44] K. Washiyama, D. Lacroix, Phys. Rev. C 78, 024610 (2008) 

[45] L. Guo, in proceeding of the NSRT12 conference (2012), to appear in Eur. Phys. J. Web of conf. 
[46] R. Bass, Phys. Rev. Lett. 39, 265 (1977) 

[47] C.R. Morton, A.C. Berriman, M. Dasgupta, D.J. Hinde, J.O. Newton, K. Hagino, I.J. Thompson, Phys. Rev. 

C 60, 044608 (1999) 
[48] A.S. Umar, V.E. Oberacker, Eur. Phys. J. A 39, 243 (2009) 

[49] F. Videbaak, R.B. Goldstein, L. Grodzins, S.G. Steadman, T.A. Belote, J.D. Garrett, Phys. Rev. C 15, 954 
(1977) 

[50] M. Evers, M. Dasgupta, D.J. Hinde, D.H. Luong, R. Rafiei, R. du Rietz, C. Simenel, Phys. Rev. C 84, 054614 
(2011) 

[51] A.S. Umar, V.E. Oberacker, J.A. Maruhn, Eur. Phys. J. A 37, 245 (2008) 
[52] C. Simenel, B. Avez, Int. J. Mod. Phys. E 17, 31 (2008) 
[53] K. Washiyama, S. Ayik, D. Lacroix, Phys. Rev. C 80, 031602 (2009) 
[54] C. Simenel, Phys. Rev. Lett. 105, 192701 (2010) 

[55] B. Yilmaz, S. Ayik, D. Lacroix, K. Washiyama, Phys. Rev. C 83, 064615 (2011) 

[56] S.E. Koonin, K.T.R. Davies, V. Maruhn-Rezwani, H. Feldmeier, S.J. Krieger, JW. Negele, Phys. Rev. C 15, 
1359 (1977) 

[57] L. Corradi, G. Pollarolo, S. Szilner, J. Phys. G 36, 113101 (2009) 

[58] R.A. Broglia, A. Winther, Heavy Ion Reactions: Lecture Notes the Elementary Processes (Frontiers in 

Physics) (Addison Wesley Publishing Company, New- York, 1991), ISBN 0201513927 
[59] W. Swiatecki, Nucl. Phys. A 376, 275 (1982) 

[60] J. Blocki, J. Randrup, W. Swiatecki, C. Tsang, Ann. Phys. 105, 427 (1977) 
[61] B. Avez, Ph.D. thesis, University of Paris XI (2009) 

[62] J. T6ke, R. Bock, G. Dai, A. Gobbi, S. Gralla, K. Hildenbrand, J. Kuzminski, W. Miiller, A. Olmi, H. Stelzer 

et al., Nucl. Phys. A 440, 327 (1985) 
[63] R. du Rietz, D.J. Hinde, M. Dasgupta, R.G. Thomas, L.R. Casques, M. Evers, N. Lobanov, A. Wakhle, 

Phys. Rev. Lett. 106, 052701 (2011) 
[64] A. Wakhle et al., in preparation 

[65] J. Reinhardt, B. Miiller, W. Greiner, Phys. Rev. A 24, 103 (1981) 

[66] E. Ackad, M. Horbatsch, Phys. Rev. A 78, 062711 (2008) 

[67] RY. Cusson, J.A. Maruhn, H. Stocker, Z. Phys. A 294, 257 (1980) 

[68] C. Simenel, C. Golabek, D.J. Kedziora, EPJ Web of Conferences 17, 09002 (2011) 

[69] V.V. Volkov, Phys. Rep. 44, 93 (1978) 

[70] C. Simenel, Phys. Rev. Lett. 106, 112502 (2011) 



